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We have studied the temperature dependence of resistivity, p, for a two-dimensional electron 
system in silicon at low electron densities, ris ~ 10^^ cm~^, near the metal/insulator transition. 
The resistivity was empirically found to scale with a single parameter. To, which approaches zero at 
some critical electron density, ric, and increases as a power To oc \ns — Wcl'' with /? = 1.6 ± 0.1 both 
in metallic (ra^ > Wc) and insulating (ra^ < Wc) regions. This dependence was found to be sample- 
independent. We have also studied the diagonal resistivity at Landau level filling factor = 3/2 
where the system is known to be in a metallic state at high magnetic field and in an insulating 
state at low magnetic field. The temperature dependencies of resistivity at 5 = and at = 3/2 
were found to be identical. These behaviors suggest a true metal/insulator transition in the two 
dimensional electron system in silicon at 5 = 0, in contrast with the well-known scaling theory. 

PACS numbers: 71.30.-|-h, 73.40.Qv, and 73.20.Fz 



I. INTRODUCTION 

For a number of years it has been generally believed 
that at zero magnetic field, all the states are localized 
in the two dimensional electron system (2DES) in the 
limit of infinite sample size. Arguments based on scal- 
ing theory indicate that as T ^ 0, resistivity always in- 
creases, exponentially in the case of "strong" localiza- 
tion or logarithmically in the case of "weak" localization 
[1]. Early experiments on relatively low mobility sam- 
ples [2,3] confirmed this behavior. Thus, conventional 
wisdom has been that there is no metal/insulator (M/I) 
phase transition in an infinite 2D sample because there 
can be no metal, in contrast with the three dimensional 
situation where electrons are localized only if the Fermi 
energy, Ep, lies below some mobility edge, E^. Recently, 
however, there have been a number of experimental and 
theoretical results which have lead us to study this prob- 
lem further. 

Experiments at high magnetic fields associated with 
the quantum Hall effect, QHE, clearly show (see, e.g., 
Ref. [4]) extended states at Ep when the Landau filling 
factor v = Ughc/eB = 1/2 (here rig is the 2D electron 
density, h is the Plank constant, c is the speed of light, 
e is the electron charge, B is magnetic field, and i is an 
integer). Recently there has been considerable interest 
in the transition from QHE behavior to insulator behav- 
ior at lower magnetic fields [5-8]. Here, theory indicates 
that the extended states will fioat up in energy as 5 ^ 
leading to an insulating state [9,10]. Experimentally, the 



case is not so clear; Shashkin et al. [7] have found recently 
that the extended states coalesce and fioat up as 5 ^ 
but do not necessarily lead to an insulator. New theo- 
retical arguments [11] have also fueled a re-examination 
of the behavior of a 2DES at 5 = 0. In Ref. [11], Azbel 
finds that for a system of noninteracting 2D electrons in 
a model disorder potential with a random set of special 
scatterers at 5 = is localized only when Ep < Eg. At 
all energies above Eg, extended states exist. According 
to Azbel, the disagreement between his results and those 
of Abrahams and coworkers [1] might indicate that the 
resistance strongly depends on the range of the scatter- 
ing centers. In the dilute 2D electron system in silicon 
metal-oxide-semiconductor field-effect transistors (MOS- 
FET's), scatterers are short-range, similar to the model 
potential used by Azbel, and thus his results may be ap- 
plicable. 

For these reasons we recently studied [12] high-mobility 
Si MOSFET samples, and showed that conventional weak 
localization, observed at T ^ 1 — 2 K, is overpowered by 
an order of magnitude drop in resistivity, p, as the tem- 
perature is decreased below ~ 1 K. No signs of electron 
localization are seen down to the lowest available temper- 
ature (20 mK) even for very low electron densities above 
some critical value, rig. At electron densities lower than 
this critical value, the resistivity monotonically increases 
as T ^ 0, indicating a localized state studied extensively 
elsewhere [13]. Moreover, the resistance is empirically 
found to scale with temperature at densities both below 
and above rig [12]. Mathematically this indicates that 
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the resistivity can be written in the form 
p{T,n,) = p[T/To{n,)]. 



(1) 



The scaling parameter, To, has been found to approach 
zero at rig = n^- The observed behaviors resemble an 
M/I transition in three dimensions rather than the be- 
havior expected for a 2D system [14], and suggest a true 
M/I transition in Si MOSFET's at zero magnetic field, in 
contrast with predictions of Ref. [1] and consistent with 
Ref. [11]. It should be mentioned that an analogous tem- 
perature behavior of resistivity and scaling on both sides 
of the transition point was observed in thin disordered 
Bi films where it was considered as evidence for a direct 
phase transition from an insulating state to a supercon- 
ducting state [15]. 

Here we report extensive studies of scaling behavior in 
high-mobility 2DES in silicon at 5 = 0. We show that 
the dependence of the scaling parameter. To, on electron 
density is symmetric about and obeys a power law 



To = C\8nf with 13 = 1.6 ± 0.1 



(2) 



where 6„ = ng —ric, both for insulating and metallic sides 
of the transition. Both the power /3 and the coefficient 
C were found to be essentially sample-independent. We 
related To for the insulating side to a localization length ^ 
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FIG. 1. Dependencies of p in units of h / 6^ on electron density for 
different temperatures at = 0. Sample Si-12b. 



using the equation ksTo = e^/e^ (where e is the dielectric 
constant) because we found p(T) cx exp (To /T) which 
is characteristic for a Coulomb gap [17]. The dependence 
of p on the electron density close to the transition point 
for different temperatures appeared to be qualitatively 
the same as those for QHE/insulator transition [6,8], 
with one temperature independent crossing point (see 
Fig. 1). We have also performed temperature studies of 
the diagonal magnetoresistivity, p^x, at = 3/2 where a 
2DES is known [4] to be in a true metallic state in the 
high ris/high B limit and in an insulating state in the 
low ris/low B limit. Therefore, at = 3/2, there must 
exist a true M/I transition observable at constant i> with 
decreasing B and as shown schematically for the phase 
diagram of Refs. [5,7,16] in Fig. 2. We have found the 
temperature and density behavior of p at 5 = to be 
absolutely identical to that for p^x at i> = 3/2, provid- 
ing additional evidence for a true metal/insulator phase 
transition in zero magnetic field. 



II. EXPERIMENTAL 

Below we show results obtained with three samples 
from two wafers: Si-12a with maximum mobility, Pmax, 
of 3.5 X 10^ cmVVs, Si-12b with prnax = 3.0 x 10^ 



cm /Vs, and Si-15 with pr. 



7.1 X 10^ cmVVs. 



Other samples from other wafers showed similar results. 
The samples are rectangular Hall bars with a source 
to drain length of 5 mm, a width of 0.8 mm, and an 
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FIG. 2. Schematically shown metal/insulator phase diagram in the 
{B J Ug^ plane taken from Refs. [5,7,16]. Lower shaded region cor- 
responds to the insulating state, two upper shaded regions — to 
the quantum Hall effect with (Txy — 1 / h and 2 6^ fh. Blank 
strips between shaded regions correspond to metallic states. The 
arrow shows the metal/insulator transition at 1^ = 3/2 realized in 
the current experiment. 
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FIG. 3. Temperature dependencies of the resistivity (sample Si- 
12b) for different electron densities at _0 = 0. 

intercontact distance of 1.25 mm. The thickness of the 
oxide separating 2DES and the gate was close to 2000 A 
for all samples. The resistance was measured using a 
four-terminal dc technique with cold amplifiers (input 
resistances > lO^'* SI) installed on the 1 K pot of a di- 
lution refrigerator. The output of these amplifiers was 
connected to a standard digital voltmeter. Great care 
was taken to ensure that all data were obtained in the re- 
gion of linear I — V characteristics. To accomplish this, 
it was necessary to chose a proper measuring current 
which varied from several pA for high resistances to 100 
nA for low resistances. This was especially important at 
low temperatures where p varies by 5 orders of magni- 
tude within the chosen interval. For each sample we 
observed the same p{T,ng) characteristics independent 
of contact configuration. Samples were mounted with a 
weak thermal link to the mixing chamber (via a stainless 
steel rod) allowing a change in the temperature from 0.2 
to 7.5 K during the experiment. We controlled the tem- 
perature using two calibrated resistance thermometers 
placed in good thermal contact with the sample. 

Figure 1 shows the resistivity (in units of h/e^) as a 
function of electron density for Si- 12b for several tem- 
peratures. One can see that all curves cross at some 
resistivity, p° ~ 2/i/e^, and electron density, = 
0.96 X 10^^ cm~^, which corresponds to a mobility of 



about 10^ cm^/Vs. At densities below this point, the re- 
sistivity is higher for lower temperatures, behavior which 
is characteristic of an insulating state. In contrast, for 
Hg > ric, the lower the temperature, the lower the re- 
sistance, behavior which is characteristic of a metallic 
state. The data, including p° , for other samples are iden- 
tical except ric varies (see also Ref. [12]). This behavior, 
particularly the existence of a single crossing point, is 
qualitatively identical to the behavior of the QHE to in- 
sulator transition (see Figs. 2 and 4 in Ref. [6] and Fig. 3 
in Ref. [8]). 

To see the temperature dependence of resistivity, in 
Fig. 3 we replot p data as a function of temperature for 
30 different electron densities varying from 7.12 x 10^*^ to 
13.7x10"'^'^ cm~^. At low densities, the curves grow mono- 
tonically as the temperature decreases, behavior charac- 
teristic of an insulator [18]. However, for ^ ric, the 
temperature behavior of p becomes nonmonotonic: resis- 
tivity increases at T > 2 K and decreases as the temper- 
ature is decreased; this behavior is "insulating" at higher 
T and "metallic" at lower T. At still higher ng, resistiv- 
ity is almost constant at T ^ 4 K but falls by an order 
of magnitude at lower temperatures showing a strongly 
metallic behavior as T ^ 0. 

A striking feature of the p(T) dependencies for differ- 
ent Hg is that they can be made to overlap by scaling 
them along the T-axis. In other words, resistivity can 
be represented as a function of T/Tq with To depending 
only on . This was possible for quite a wide range 
of electron densities (typically — 2.5 x 10"'^'^ ^ '^s ^ 
ric + 2.5 X 10"'^'^ cm~^) and in the temperature interval 
0.2 to 3 K. The results of this scaling are shown in Fig. 4 
where p is represented as a function of T/Tq. One can 
see that the data dramatically collapse into two sepa- 
rate curves, the upper one with open symbols for the 
insulating side of the transition and the lower one with 
filled symbols for the metallic side. The thickness of the 
lines is largely governed by the noise within a given data 
set which indicates the high quality of the scaling. It is 
worth noting that qualitatively the same scaling picture 
was obtained in Ref. [15] for insulator/superconductor 
transition in thin Bi films. 

The procedure used to bring about the collapse and 
determine To for each was the following. First, we 
noted (see Fig. 5) that at low n^, to high precision the 
resistivity has a temperature behavior characteristic of 
hopping in the presence of a Coulomb gap [17]: 

p = Poexp[(To/T)i/'] . (3) 

Therefore, Tq^^ for the lowest (i.e., for the upper curve 
in Fig. 3) was determined by fitting Eq. 3 to the resis- 
tance data. Next, the second curve from the top in Fig. 3 
was scaled along the T-axis with the factor 7^^) to coin- 
cide with the upper curve; To for the second curve thus 
was determined as Tg^-* = j''-^'^Tq^\ Then the third curve 
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FIG. 4. Resistivity vs T/Tq (a) and scaling parameter, Tq, vs 
electron density (b) for Si-12a. Open symbols correspond to the 
insulating side of the transition, closed — to the metallic one. 

from the top was scaled with the factor 7^^) to coincide 
with the two upper ones determining Tg^-* = j''-^'^Tq^\ 
This procedure was repeated subsequently for all "in- 
sulating" curves (rig < ric, dp/dT < 0) in Fig. 3. This 
gives the upper curve in Fig. 4 (a) designated with open 
symbols. The scaling parameter, To, for the insulating 
side of the transition, is shown in Fig. 4 (b) by open sym- 
bols as a function of electron density. One can see that 
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FIG. 5. Resistivity vs T for two samples at lowest electron 

densities which were achieved. 

it approaches zero as approaches the critical point, 
rig. Identical behavior for a second sample is shown in 
Fig. 6. 

To collapse the "metallic" curves (rig > rig) into a sin- 
gle curve, we applied the same procedure. We started 
with the curve corresponding to the highest (i.e., with 
the lowest curve in Fig. 3) and then scaled the T-axis for 
the other curves by factors Y*-^-*, Y*-^-* etc to make them 
coincide with the first curve. To fix the position of the 
resulting curve in the log(T/To) axis, it is necessary to 
assign some value of To to the first curve, and then get 



To for other curves as T 



(k) 



1 



(*^)Tn^^\ In contrast to 



the "insulating side" , it is not clear what value of To 
should be assigned to the first metallic curve. We have 
found (see Fig. 7) that for both the insulating side and 
the metallic side of the transition, the dependence of 7 
on \bn \ = \ns — ricl is a power law: 



(4) 



with the same value of [3 for metallic and insulating 
sides. This common power law can be clearly seen in 
Fig. 7 where for each sample the open (insulating side) 
and filled (metallic side) symbols form a single line. Be- 
cause 7(|i5n|) has the same behavior on both sides of the 
transition, we suggest that the function To(|(5„|) is also 
quantitatively symmetric. In Fig. 10 and in the discus- 
sion below we further justify this assumption [19]. The 
scaling factor as a function of electron density for the 
metallic side of the transition is shown in Figs. 4 (b) and 
6 (b) by closed symbols. Again, as for the insulating side. 

To ^ as ris ^ ric- 

Figure 7 shows To (in log-log format) for both metallic 
and insulating sides of the transition for three different 
samples as a function As we mentioned above, the 

dependency To(|(5„|) is a power law. We find the aver- 
age power [3 to be 1.60 ± 0.1 for insulating side of the 
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FIG. 6. The same as Fig. 4 but for Si-12b. 

transition and 1.62 ± 0.1 for metallic side. For all three 
samples shown in the figure, the dependencies To(|(5„|) 
are nearly identical in spite of the fact that the values of 
ric and sample mobilities are different. We should men- 
tion that the similar power law with an exponent 1.5±0.2 
was observed for scaling of the superconductor/insulator 
transition in thin Bi films [15,20] though we believe the 
physical mechanism driving the transition is different in 
the two systems. 

Finally, we show the temperature dependencies of the 



10^ 



10' 



10' r 



10" 



10"' r 



10"' 





circles - Si-15 




squares - Si-12b 


• 

• 


diamonds - Si-1 2a 








y , * 




# .* 
























♦ 

• 




□ 




h ■ 









10 10 

I 5 I (lO" cm"') 



10" 



10 



FIG. 7. Scaling parameter, Tq, as a function of \ for three 
samples from different wafers. 
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FIG. 8. Resistivity vs temperature for Si-12a in a magnetic field 
corresponding to = 3/2 (symbols) and at _0 = (dashed 
curves) for different electron densities. 

diagonal resistivity, p^x, in a magnetic field correspond- 
ing to a Landau level filling factor i> = 3/2. To obtain 
this data, we varied both and B so that i> remained 
constant as shown in Fig. 2. Along this path one expects 
a true M/I transition. Figure 8 shows Pxx(T) dependen- 
cies for i> = 3/2 (symbols) along with p(T) dependencies 
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FIG. 9. Resistivity vs T/Tq for Si-12a. Open symbols — the 
insulating side of the transition; closed — the metallic one. 

at 5 = (dashed lines). Surprisingly, the curves are 
practically identical even though is different. In par- 
ticular, curves with high-T resistivities of similar values 
were found to have identical temperature dependencies 
indicating that p'-' is also the same. This shows that 
the M/I transitions at zero magnetic field and at half- 
integer filling factor are identical when viewed in terms 
of resistivities. 



III. DISCUSSION 

Scaling of an appropriate physical variable is one of the 
hallmarks of a phase transition. In this paper we report 
the scaling behavior for the resistivity of the high-quality 
2DES in Si MOSFET's with the scaling parameter ap- 
proaching zero as the density of electrons approaches the 
critical value. This behavior strongly suggests a true 
metal/insulator transition in the two-dimensional sys- 
tem in silicon at zero magnetic field, in apparent con- 
tradiction to the theoretical arguments of Abrahams and 
coworkers [1] and consistent with those of Azbel [11]. 
We must note however that the temperature behavior 
of resistivity on the metallic side of the transition with 
p dropping by an order of magnitude at temperatures 
below 1-2 K is not what one would expect for an or- 



dinary metal where resistivity saturates when the fre- 
quency of phonons = ksT/h becomes less than r~^, 
the inverse elastic scattering time, i.e., at temperatures 
T "^Tq = h/ksT (in the regime of interest, Tq > 10 K). 
We have discussed in detail the precipitous drop of re- 
sistivity at T < 1 — 2 K in Ref. [12]. In particu- 
lar, we showed that this drop can be explained neither 
by electron-phonon scattering nor by the temperature- 
dependent screening considered in Refs. [21,22]. 

In Fig. 9 we plot the resistivity as a function 
of [To/T(ng)Y^^ . The resulting "metallic" curve is 
shown by closed symbols (lower curve). It mono- 
tonically decreases by more than an order of magni- 
tude as [To/T(ng)Y^^ increases and finally saturates at 
[To/T(ng)Y^^ ^ 6. The upper, "insulating" curve (open 
symbols) at (T/To)-^/^ > 2 can be fitted well by the 
formula 



p = Aexp[(T/To)-i/2] 



(5) 



where the prefactor A is close to h/e^. At (Tq/TY^^ ^ 2 
(i.e., when T becomes close to To) the dependence weak- 
ens. As we have already mentioned, this behavior is 
characteristic of hopping in the presence of a Coulomb 
gap [17]. This contradicts earlier experimental results 
of Ref. [23] where the absence of the Coulomb gap in 
Si MOSFET's was reported. We attribute the differ- 
ence to the greater role of Coulomb effects relative to 
scattering by impurities in our samples because our sam- 
ples are of higher mobility (and therefore are less dis- 
ordered) at low electron densities. Recently, Coulomb 
gap behavior of the resistivity was reported for low- 
density gated GaAs/AlGaAs heterostructures [24]. In 
this system, a transition from the Coulomb gap behav- 
ior to the Mott variable range hopping behavior [25], 
p (X exp [(T/To)~^^^], was observed as the temperature 
was lowered, in agreement with recent theoretical predic- 
tions [26]. In principle, in our system, p(T) characteristic 
of the Mott variable range hopping regime can develop at 
lower temperatures and lower electron densities leading 
to the destruction of scaling as T ^ 0. 

Qualitatively the same p(T/To) function was observed 
for the insulating state in thin palladium films [27]. The 
resistivity for the thinnest films also was found to obey 
the Coulomb gap law. However, the system studied in 
Ref. [27] was always weakly or strongly insulating; there 
was no "metallic part" of the scaling function. 

We obtain physical insight into the observed scaling 
by relating the scaling parameter. To, to a length scale 
as indicated in Ref. [17], 



(6) 



Therefore, if we assume that e is constant 
[28], the approach of To to zero as 6„ ^ 
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FIG. 10. Resistivity vs the scaling variable, jT^^ '^^ with Z — \ 
and = 1.6. Open symbols — the insulating side of the transi- 
tion; closed — the metallic one. 



indicates a diverging localization length: 
^ (X with P ^ 1.6 



(7) 



(the conventional symbol for this exponent is i>, but we 
have used /3 in this paper to avoid confusion with i>, the 
Landau filling factor). Similar behavior for the local- 
ization length (with /3 1) was recently deduced from 
nonlinear current- voltage characteristics [16]. 

Figure 10 shows p as a function of the scaling vari- 
able, X = bn/T^I'P, used in Ref. [29] for the descrip- 
tion of the superconductor/insulator transition. Here z 
is the dynamical exponent equal to 1 in the case of purely 
Coulombic interaction with constant e. The exponent [3 
was taken to be 1.6. Again, all the data collapse onto 
a single curve. One can see that p as a function of the 
scaling variable, X , is continuous and smooth around the 
transition point, X = Q (see inset). This indicates that 
the mechanism responsible for the temperature depen- 
dence of resistivity on both insulating and metallic sides 
of the transition is the same, and justifies our assumption 
about qualitative symmetry of the scaling parameter. To, 
around the critical point. 

Let us now discuss Pxx(T) for filling factor i> = 3/2. 
In contrast with the zero magnetic field situation, there 
is no doubt that in high enough magnetic fields there are 
extended states at the center of each Landau level [4] 
corresponding to Landau level filling factors = i -|- 1/2. 
At low magnetic fields, the extended states no longer fol- 
low the centers of the Landau levels but are expected to 
fioat up in energy [9,10]. It was recently shown experi- 
mentally [7] that the extended states indeed fioat up in 
energy at low magnetic fields. In contrast with theoreti- 
cal expectation, however, their energy does not increase 
infinitely; instead they combine at some finite energy. An 
insulator/metal/QHE phase diagram for high-mobility Si 



MOSFET's (similar to those used in this work) was ob- 
tained in Refs. [5,7,16], and is shown schematically in 
Fig. 2. One can see metallic strips at = 1/2 and 3/2, 
corresponding to the extended states at these filling fac- 
tors, between quantum Hall states with a^y = le^/h 
and 2e^/h and insulating state with axy = (metallic 
strips at other half-integer filling factors are not shown). 
Staying at the same filling factor i> = 3/2, therefore, 
one can observe a transition from metallic state in the 
high 5/high limit to the insulating state at lower B 
and Hg. 

Temperature dependencies of the diagonal resistivities 
shown in Fig. 8 for both i> = 3/2 (symbols) and zero 
magnetic field (dashed lines) are essentially identical. 
The curves having close "starting resistivities" at high 
temperatures have the same temperature dependencies. 
This strongly suggests that the M/I transition in zero 
magnetic field is identical to the M/I transition at half- 
integer filling factor. We note that the identical M/I 
transitions at 5 = and at high B were indicated in 
Ref. [16] where nearly field-independent behavior of the 
localization length was reported. 



IV. SUMMARY 



In summary, we have shown that the zero magnetic 
field resistivity of the high-quality two-dimensional elec- 
tron system in silicon scales with temperature. A single 
scaling parameter collapses the resistivity data onto two 
curves, insulating for < rig and metallic for > ric, 
and decreases upon approaching the critical electron den- 
sity as a power /3 = 1.6±0.1 for both metallic and insulat- 
ing sides of the transition. This scaling behavior strongly 
suggests a metal/insulator phase transition for the 2DES 
in Si MOSFET's at 5 = 0. We have compared this tran- 
sition with the phase transition in a magnetic field corre- 
sponding to a constant Landau level filling factor i> = 3/2 
and with the superconductor/insulator phase transition 
in thin Bi films. The behavior we observe is identical 
to the i> = 3/2 case and quite similar to the Bi case. 
These similarities add further weight to our identification 
of the reported transition as a true M/I phase transition. 
In particular, the identical behavior at = 3/2 and at 
5 = along with the arguments presented in Ref. [12] 
indicate that the observed scaling behavior is not due to 
some sample size effect associated with weak localization. 
Thus, our data appears to be consistent with the work 
of Azbel [11], while the conclusions of Ref . [1] seem inap- 
propriate for high-mobility Si MOSFET's. We are in the 
process of extending our measurements to GaAs/AlGaAs 
samples, of studying the scaling behavior at low B and 
with electric fields. 
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